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Abstract 

An alternative approach to the Jaynes-Cummings model (JCM) with dissipation 
at a finite enviromental temperature is presented in terms of a new master equation 
under Born-Markovian approximations. An analytic solution of the dissipation JCM 
is obtained. A variety of physical quantities of interest are calculated analytically. 
Dynamical properties of the atom and the field are investigated in some detail. It 
is shown that both cavity damping and environmental temperature strongly affect 
nonclassical effects in the JCM, such as collapse and revivals of the atomic inversion, 
oscillations of the photon-number distribution, quadrature squeezing of the field and 
sub-Poissonian photon statistics. 



1 Introduction 



In the past few years there has been considerable interest in studying the Jaynes-Cummings 
(JCM) [1] with dissipation since the experimental verification of the collapse and revivals 
of Rabi oscillations [2-4]. A number of authors have treated the JCM with dissipation 
by means of analytic approximation [5,6] as well as numerical calculations [7-11]. Most 
of these work dealt with the JCM with dissipation only at the zero enviromental tem- 
perature due to some technical difficulties. As is well known, in the JCM one must solve 
a master equation with cavity damping. Generally speaking, it is more difficult to solve 
the master eqution analytically. To our knowledge, no analytic solution for this model 
has been given for general initial conditions and environmental temperatures. At zero 
temperature, Agaral and Puri [12] presented an analytic solution for the initial state of 
the field being a vacuum state; And Daeubler et al. [13] found an analytic expression 
for the atomic inversion and the intensity of the cavity field for the initial state of the 
field being a coherent state by using the quasiprobability-distribution technique [14]. At 
a finite temperature, only some numerical solutions [8,9,11] have been reported. The 
purpose of the this paper is intend to present an alternative approach to the dissipation 
JCM at a finite temperature by introducing a new master equation which comes from a 
simple coupling between the system and its environment. We study dynamical properties 
of the field in the JCM analytically. In particular, we show that both cavity damping and 
environmental temperature strongly affect nonclassical effects in the JCM, such as col- 
lapes and revivals of the atomic inversion, oscillations of the photon-number distribution, 
quadrature squeezing of the field, and sub-Poissonian photon statistics. 

This paper is organized as follows: In Sec. 2, a new master equation is presented and an 
explicit expression of the master equation for the JCM is obtained. In Sec. 3, dynamical 
properties of the atom is investigated. The atomic inversion, the dipole moment of the 
atom and the atomic entropy are calculated analytically. Sec.4 is devoted to dynamics 
of the field. Oscillations of the photon-number distribution, quadrature squeezing of the 
field and sub-Poissonian photon statistics are discussed in detail. Finally, in Sec. 5 we will 
summarize our results and give some concluding remarks. 
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2 Master Equation and its Analytic Solution 
for the JCM 



Let H and p to be the Hamiltonian and the density operator of the JCM, respectively. 
We use a bath of harmonic oscillators to model the cavity damping which describes the 
irreversible motion of p caused by the enviroment consisting of an infinite set of harmonic 
oscillators. We assume that the atom-field system of the JCM the system interacting with 
the environmemt (or a reservoir) can be described by the total Hamiltonian 

Ht^H + Y. ^^Mh + ^HY: C{u;,) ik + H) + h^H^ E (2-1) 

where hi and 6^ are the boson annihilation and ceration operators for the enviromemt, in 
Eq.(2.1) the second term is the Hamiltonian of the reservoir, the third one represents the 
interaction between the system and the reservoir , and the last one is the renormalization 
term which compensates for the coupling-induced renormalization of the potential [15,16]. 
In the total Hamiltonian we have adopted a simple coupling between the system and the 
reservoir such that it satisfies the condition \y,H\ — {){V denotes the interaction term in 
(2.1)), which is required in the back-action evading and quantum-nondemolition schemes 
[17] and applied to decoherence of quantum system [18-20]. 

Following the standard master-equation approach [21], making use of the Hamiltonian 
(2.1) we can derive the following master equation for the reduced density operator in the 
Schrodinger picture under Born-Markovian approximations: 

^ = ^[H, pit)] - l[H, [H, Km - ^^[H, Pit)] (2.2) 

where 

, , kT ^. Jiu)\Ciu)\^ 
7 = Au;' + — lim ^ ^' ^ ^' , (2.3 

Jiu;)\Ciuj 



^0 UJ 



with 



^ , , , J(w)|C(a;)P 

Ao;' = ih / do; ^ ^' ^ ^' (2.5) 



10 UJ 

In the above equations, Jiui) is the spectral density of the reservoir, V is the Cauchy 
principal part of the integration [21], T is the temperature of the enviroment and k is the 



Boltzmann constant. In the derivation of the master equation, we have assumed that the 
temperature T is high enough so that the Markovian approximation is vahd. 

If we neglect the the Lamb shift term, i.e., the last term on rhs of Eq.(2.1), the master 
equation becomes 

^ = l[^,p(t)]-7[^,[^,pW]] (2.6) 

Notoce that this equation has the same form as the Milburn's equation [22] under diffu- 
sion approximatiom. However, they come from completely different physical mechanism. 
The former originates from the dissipation while the latter is from the uncontinuous and 
stochastic unitary evolution. In what follows we will use the master equation (2.6) to 
study the JCM with dissipation. 

The resonant Jaynes-Cummings Hammiltonian describing an interaction of a two-level 
atom with a single-mode cavity field in the rotating-wave approximation is given by 

H = huj{a^a + —0-3) + %\{a_a^ + a+a) , (2.7) 

where cu is the resonant frequency of the cavity field and the atomic trasition, A is the 
atom-field coupling constant; o and o"*" are the field annihilation and creation operators, 
respectively; (T3 is the atomic-inversion operator and a± arc the atomic "spin flip" oper- 
ators which satisfy the relations [(T+,a"_] = 2(T3 and [(3"3,a"±] = ±2a±. For simplicity, we 
take h — 1 throughout the paper. 

We now find the analytic solution of the master equation (2.6) apphed to the Hamil- 
tonian (2.7). Moya-Cessa et al. [23] have obtained an formal solution of (2.6) for a 
Hamiltonian with a little diffrence from (2.7). In what follows we will present an explicit 
solution of the master equation for the Hamiltonian under our consideration. Following 
the approach in refs. [23,24], we introduce three superoperators R, S and T which are 
defined through their actions on the density operator, respectively, 

exp{Rt)m = J^^-jT^H'^mH' (2.8) 

exp(5t)p(0) = exp(-i^t)p(0) exp(i^t) (2.9) 
exp{ft)p{Q) = exp(-7t^2)p(0) exp(-7t^2) (2.10) 

where p(0) is the initial density operator of the atom-field system, and the Hamiltonia H 
is given by Eq.(2.7). 
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It can be checked that the master equation (2.6) has the following formal solution: 

p{t) = exp(^t) exp(5f) exp(f t)p(O) (2.11) 
We assume that initially the field is prepared in the coherent state | z) defined by 

oo 1 

\z) = Y.Qn\n), Qr, = exp{--\z\')^ (2.12) 

n=0 ^ V*^- 

and the atom was prepared in its excited state | e), so that the initial density operator of 
the atom-field system takes this form: 

m = ('"^^o^"' o) (2.13) 

We divide the Hamiltonian (2.7) into a sum of two terms which commute with each 
other, i.e., 

H^Ho + Hi, [Ho,Hi]^0 (2.14) 

with 

where n — a^a. 

Similarly, we can make the following decomposition: 

H^^A + B, [A,B]^0 (2.16) 

where the representations of the operators A and B in the two-dimensional atomic basis 
take the following forms: 

i = f M B = 2Xu;( V.^ - ^) ] (2.17) 

where 

5-^ = cu^n - + X^h (2.18) 
For convenience, we introduce the following auxiliary operator: 

P2{t) = exp{St) exp(f t)p(O) (2.19) 
Prom the definition of the superoperators and the initial condition (2.13), we find that 

P2{t) = ex.p{-iHit) ex.p{-'ytB)pi{t) exp(-7i5) ex.p{iHit) (2.20) 



Here the operator pi{t) is defined by 
^(0) ^ ( I W I I ) (2.21) 



where 



I ^(t)) = exp{-^t[u\h + + A2aa+]} | ze''^') (2.22) 
For the exponential operators on the rhs of Eq.(2.20) we can find that 



exp(-z^,t) = . (2.24) 



where 

^^(t) = cos(Atv^), Snit) = sm{XtVh) (2.25) 

Xn{t) = cosh[2A7a;t(ri - ^)v^], Fn(*) = sinh[2A7a;t(n - ^)v^] (2.26) 

Then, from Eqs.(2.23) and (2.24) it follows that 



(2.27) 



where 

Rn = Cn{t)Xn{t) + tSn{t)Yn{t) (2.28) 

K = Cnit)Yr,{t) + iSn{t)Xr,it) (2.29) 

Substituting Eq.(2.27) into Eq.(2.20), we can obtain an explicit expression for the 
operator p2{t) as follows: 

where we have used the following symbol: 

%{t)^\%{t)){^j{t)l (i,j = l,2) (2.31) 

with 

I ^i(t)) = i?„+i(t)|*(t)), I ^,{t) = -M)a+|vl/(t)) (2.32) 
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where | "^{i}) is given by equation (2.22). 

Taking into account the definition of the superoperator R, through the action of the 
operator exp{Rt) on the operator P2{t) one can obtain the following formal solution: 

Pit) = f:^^H'p,{t)H' (2.33) 

fc=0 

where the operator H and P2{t) are given by Eqs.(2.7) and (2.30), respectively. 

Indeed, Eq.(2.33) is the exact solution of the master equation (2.6) for the resonant 
Jaynes-Cummings Hamiltonian (2.7) in the operator form. Although the form of the 
solution (2.33) is pleasant, it is unconvenient in use. In most cases of practical interest, 
one needs to know the matrix elements of the density operator p{t) in the two-dimensional 
atomic basis to calculate expectation values of observablcs. 

In order to find the explicit form of the solution, we need that 



where 



fn{k) = ^(<^^ + fy, Uk) & (2-35) 



with 



1 - 1 

ipn = (^(^ ~ 7^) + -^V^, (j)n = i^in — -) — A\/^ (2.36) 

For convenience, we define the following matrix: 

M^^\t)^H^P2{t)H^ (2.37) 

Prom Eqs.(30) and (2.34) we can obtain its matrix elements as follows: 

Mfiit) = /n+i(A:)^ii(t)/n+i(/c)+a^;(A;)^2i(t)/n+i(fc) 

+/n+i(A;)*i2(i)^;(/^)a+ + aUk)^22{t)UkYa^ (2-38) 
M^Sit) = ^;(A;)o+^n(t)a^;(A;) + /„(A;)^2i(t)a^;(A;) 

+g'Ma+^i2{t)fn{k) + fn{k)'^22{t)fn{k) (2.39) 

- g'^{k)a^^^^{t)dU+^{k) + fn{k)4'2i{t)U+i{k) 

+g'M^^^i2{tmk)a^ + Uk)i>22{t)Ukya^ (2-40) 
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where 

tik) = -^9n{k) (2.41) 

Therefore, the analytic solution of the master equation (2.6) for the Jaynes- Cummings 
Hamiltonian (2.7) can be expressed explicitly as follows: 



Making use of this solution, one can evaluate mean values of operators of interest. In 
the next two sections we will use it to investigate various dynamical properties of the 
dissipation JCM. 

3 Dynamical Properties of the Atom 

In this section, we shall study dynamical properties of the atom and discuss the influence 
of the dissipation on them. We will derive analytic expressions of the atomic inversion 
and the dipole momentum of the atom, and calculate the atomic entropy. 

3.1 Collapse and revivals of the atomic inversion 

The atomic inversion is defined as the probability of the atom being the excited state 
minus the probability of being the ground state, that is 

W{t) = Tr[pA{t)&s] (3.1) 

where pAit) is the reduced density operator of the atom, it can be obtained through 
tracing over the field part in (2.42). Making use of the solution (2.42), one can rewrite 
the inversion (3.1) as 

= E ^[(^I-MSiV) - {n\M^M] (3.2) 

k,n=0 

We now calculate the two expectation values on the rhs of the above equation. Making 
use of Eqs.(2.38), (2.39) and (2.42) we find that 

{n\Mfht)\n) = {U,{k)f\i^Mt)?H9n+i{k)f\Un + l.t)\' 

+2Re{fn+i(k)gn+i(k)ri(n,t)Mn + ht)} (3.3) 

{n\M^S{t)\n) = {gMY\Mn^ht)\'+{fn{k)Y\Mn,t)\' 

+2Re{U{k)gn{k)Mn,t)ri{n - l,t)} (3.4) 



where functions fn{k) and gn{k) are given through replacing the number operator fi by 
the number n in Eq.(2.35), and we have introduced the following symbols: 

iPi{n, t) = {n\ ^i(t)), ^2(n, t) = {n \ ^2{t)) (3.5) 

which can be explicitly written as 

= ^Qn |c„+i[l + exp{A-fXu;t{n + ^)Vn + 1)] 

-iSn+i[l - exp{A^Xijt{n + ■ exp{-jtipl^,)e-''^' (3.6) 

ip2{n, t) = ^Qn |-C„[1 - exp(47Aa;t(n - ^)Vri)] 

+iSn[l + exp(47Aa;t(n - ^)V^)]| • exp(-7t(^2^e-*("-i)"* (3.7) 

where C„(t), Sn{t), and are obtained from their corresponding operator form through 
the replacement: h ^ n. 

From Eqs.(3.3), (3.4), (3.6), and (3.7) it follows that 

{n\Mf^{t)\n) = i|g„|2{^2fc^iexp[-27t^^+i] + </'f+iexp[-27t(/*^^i] 
+2(^^+10^ cos[2AiVn + 1] exp[-27t(^„+i0„+i] 

•exp[-4AV(^ + l)]} (3.8) 

(n|^g)(t)|n) = ^|gn-ir{^f exp[-27tv^^] + 0f exp[-27t(/<^] 
-2<^^+i(^^cos[2Atv^ exp[-27i^„0„] 

•exp[-4A27tri]} (3.9) 

Substituting Eqs.(3.8) and (3.9) into Eq.(3.2), after summing over k, we arrive at the 
result: 

oo 

W{t) = 5^ I Q„ |2 exp[-4A27t(n + 1)] cos(2AiVriTT) (3.10) 

n=0 

which indicates that the revivals of the atomic inversion are destroyed in the time evolution 
due to the appearance of the decay factor exp[— 47A^t(n + 1)] which comes from the 
contribution of the damping term in the master equation (2.6). Obviously, when 7=0, 
the usual result can be recovered. 
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3.2 The dipole moment of the atom 

As is known, the dynamical properties of the atom involve not only the atomic inversion 
but also the knowledge of the dynamics of coherence between the two atomic levels, which 
can be described by the dipole moment defined by 

D{t) = Tr[pA(t)a4 (3.11) 

which can be expressed explicitly as follows: 

m = E ^^{n\M^S\n) (3.12) 

k,n=0 

It is not trival to evaluate the mean value on the rhs of the above equation. From 
Eq.(2.40) it follows that 

{n\M^S{t)\n) = fn+i{k)fn{k)r2{n,t)Mn,t) 

+fn+i{k)gn{k)tp*i{n- l,t)tpi{n,t) 
+gn{k)gn+i{k)i}l{n - l,t)ijj2{n + l,t) 

+U{k)gn+i{k)r2{n,t)il;2{n + l,t) (3.13) 
Taking into account Eq.(2.35), one can express (3.13) as a more useful form: 

{n\M^^2\t)\n) = ^(^^+i^^[V't(n-l,i)V'i(n,i) + ^;(n,t)V^i(n,t) 

+V^*(n - 1, t)'ilj2{n + l,t)+ '0;(n, t)V'2(ri + 1, t)] 
1 
4 

-ipl{n - 1, t)tp2{n + l,t)+ i/j^in, t)ij2{n + 1, t)] 

-ipl{n - 1, t)ilj2{n + l,t)- tp^in, t)ij2{n + 1, t)] 

+\(t>l+M-i^l{n - 1, t)Vi(n, t) + V;(n, t)V'i(n, t) 

+ri{n - 1, t)V'2(n + 1, t) - V2*(^, t)i^2{n + 1, t)\ (3.14) 

Substituting Eq.(3.14) into (3.12) and making use of Eqs.(3.6) and (3.7), after summing 
over k and through a lengthy calculation we obtain the result: 



+ -^n+M-Vi{n - 1, t)V'i(n, t) + ^2*(n, t)^i(n, t) 



D{t) = -EQnQ;-iexp{-7t[a;-Aa_(n)]2} 
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exp[47Acc;t(na+(n) - -a_(n))] • e-*t^''-(")+^l* 

2 



-exp[-4jXujt{{n+^)a_{n) + \^n{n+l)] ■ e'l^'^+W"^!* 
+ exp[47Aa;t((n — -)a_(n) — 

- exp[-47Aa;t(ria+(n) + ^a_(n))] • e-^[^«-(")-^l*| (3.15) 

where 

a±(n) = x/n ± Vn + 1 (3.16) 
3.3 The entropy of the atom 

It is well known that if we assume that initially the atom and the field in the JCM are in 
a pure state,then at f > 0, the atom-field system evolves into an entangled state. In this 
entangled state the field and the atom separately are mixed states. Phoenix and Knight 
[25] have shown that entropy is the most appropriate quantity of measuring the purity of 
the quantum state in the JCM. The time evolution of the atomic (field) entropy reflects the 
time evolution of the degree of entanglement between the atom and the field. The higher 
the entropy is, the greater the entanglement between the atom and the field becomes. 

The atomic entropy is defined in terms of the reduced density operator of the atom in 
this form: 

SA{t)^-Tr{pA{t)\^PA{t)\ (3.17) 



where the reduced density operator of the atom can be written as 

All A 
A21 A22 



PA{t) = ^ (3.18) 



where 

>^^^ = E ^^{n\Mf\n), (i,J = l,2) (3.19) 

fe,n=0 

It is easy to see that the off-diagonal element A12 equals the diploe moment of the atom, 

i.e., A12 = D{t). Making use of Eqs.(3.8) and (3.9), from Eq.(3.19) one can find that 

00 

2 



All = ^ ^ I Q„ r {1 + exp[-4A27i(n + 1)] cos(2AtVn + 1)} (3.20) 



n=0 



-1 00 

A22 = - E I Qn-i r {1 - exp[-4A27H cos(2Atv^)} (3.21) 

Zi 



n=0 
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Since the trace is invariant under a similarity transformation, one can go to a basis 
in which the reduced density operator of the atom is diagonal. Then, the atomic entropy 
(3.17) can be expressed as follows: 

S'^(t) = — q;_i_ In q;_|_ — q;_ lnQ;_ (3.22) 

where 

«± = ^{1 ± ^l-4(AnA22- |Ai2|2)} (3.23) 

4 Dynamical Properties of the Field 

As is known, the field in the JCM can exhibit a lot of nonclassical effects in the time 
evolution, such as oscillations of the photon-number distribution, quadrature squeezing 
and sub-Poissonian photon statistics. In this section, we shall investigate dynamical 
properties of the cavity field in the presence of dissipation, and discuss the influence of 
the dssipation on these nonclassical effects. 

4.1 Oscillations of the phot on- number distribution 

It is easy to show that the probability p{n, t) of finding n photons in the cavity field is 
given by 

pM = E ^[(^l-A^SiV) + {n\M^M] (4.1) 

k=Q 

where two expectation values on the rhs of the above equation have been given explicitly 
in the previous section. 

Substituting Eqs.(3.8) and (3.9) into (4.1), after summing over k we obtain 

p(n, t) = ^ I <5n r {1 + exp[-4A27i(n + 1)] cos(2AiVn + 1)} 

+^ I Qn-i 1^ {1 - exp[-4A^7in] cos(2Atv^)} (4.2) 

With the help of this probability distribution, it is straightforward to obtain the mean 
number of photons in the cavity field with the result: 

1 1 - °° , 

(n(t)) = n + - - -e"" ^ — exp[-4A27t(n + 1)] cos(2AtVriTT) (4.3) 

2 2 n. 

where n = |zp is the initial mean photon number in the field. As expected, from ex- 
pressions (4.2) and (4.3) we can see that there is two decay factors expf— 4A^7t(n -|- 1)] 
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and exp[—4X^'jtn] which come from the damping term in the master equation (2.6). The 
oscillatory behaviors of p{n,t) and {n{t)) are weakened with increasing of the parameter 
7. Since the damping is proportional to the temperature of the environment, these oscil- 
latory behaviors deteriorate with the increasing of the environmental temperature. When 
the damping vanishes, (4.2) and (4.3) reduce to 



p{n,t) = \Qn\^cos^{XWn + l) + |Q„_i|^sin^(Atv^) (4.4) 
{n{t)) -n+-- -e-^ ^ — cos(2AiV'^) (4.5) 

^ ^ n=0 

which are just the usual expressions without dissipation. 
4.2 Squeezing properties of the cavity field 

We now study the quadrature squeezing of the field in the cavity field. We introduce the 
two slowly varying Hermitian quadrature components of the field and X2 defined by, 
respectively, 

= l(ae'"* + a+e"'"'*), X2 = ^(ae'"* - a+e"^"*) (4.6) 

The commutation of Xi and X2 is [^1, ^2] = f. The variances {(AXif) = {Xf) - ((X^))^ 
(i = 1,2) satisfy the Heisenberg uncertainty relation ((AXi)^) ((AXi)^) > ^. A state 
of the field is said to be squeezed when one of the quadrature components Xi and X2 
satisfies the uncertainty relation ((AXj)^) < |. The degree of squeezing can be measured 
by the squeezing parameters [26] Si {i — 1,2) defined by 

^ _ {{AX,)')-l\{[X„X2])\ 

il([^i,^2])i 

which can be expressed in terms of the annihilation and creation operators of the field as 
follows: 

Si = 2(a+a) +2Re(a2e(*'"*)) -4(Re(ae*"*))2 (4.8) 
5-2 = 2{a+a) -2Iie{a^e^'^^'^) -4{lni{ae'^')f (4.9) 

Then, the condition for squeezing in the quadrature component can simply be written as 
Si<Q. 



12 



In principle, expectation values for any function F{a'^, a) are calculated by the follow- 
ing formula: 

OO OO o 

{F{a+,a)) = Y.EiT^'k\[{n\M^i;>{t)F{a\a)^^^^ (4.10) 

n=0 ifc=0 ^ 

However, it is generally not an easy matter to calculate expectation value of an arbitrary 
function F(o+,a) for the field in the JCM. For the expectation value (ae*'^*), through a 
tedious calcutation we find that 

-1 OO 

* n=0 

• {a+(n)Ti exp{2-ftipn+ifn) + a-{n)T2exp{2jtipn+i(j)n) 

+a+{n)T3 ex.p{2-ft(fin4)n+i) + a+{n)T4 exp(27t0„+i0„)} (4-11) 

where S'„, (fn, and 0„ are given by the replacement: n — >• n in Eqs.(2.18) and (2.36), and 

Tl = [Cn+l{t)-iSn+imCnit)+iSnit)] 

■[X^+i{t) - iY^+,mXr,{t) - lYr,{t)] (4.12) 

T2 = [Cn+l{t)-tSn+l{t)][Cn{t)-zSn{t)] 

■[Xn+lit) - iYn+imXnit) + lYnit)] (4.13) 

Ta = [Cn+lit)+iSn+imCnit)+iSnit)] 

■[Xn+l{t) + tYn+imXn{t) - tYn{t)] (4.14) 

n = [C^+l{t)-lSr,+,mCn{t)-tSn{t)] 

•[X„+i(t) + iF„+i(t)][X„(t) + iY„{t)] (4.15) 

where the functions C„(t), Sn{t), Xn{t), and l^(t) are given by the replacement:??. — >• n in 
Eqs.(2.25) and (2.26). 

Substituting the exphcit expressions of the functions C„(i), Sn{t), Xn{t), and F„(t) 
into Eqs.(4.12) ~ (4.15), from Eq.(4.11) we find that 

2 CO 

{ae'^^) = -^Q„(5;_i{a+(n)exp[iAa_(n)t]exp[-6_(n)7i] 

n=0 

-|-a_(n) exp[— zAa_(n)i] exp[— 6+(n)7t] 

-|-a_(n) exp[iAa+(n)t] exp[— c_(n)7t] 

-|-a+(n) exp[— zAa+(n)i] exp[— c+(n)7i]} (4-16) 
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where 

b±{n) = [co± Aa_(r^)]^ c±(n) = [to ± Xa+{n)f (4.17) 

Similarly, we can find that 

1 CXD 

n=0 

•{a+(n)7;'exp(27t(/7„+i(^„_i) + a'_{n)T2exp{2jt(pn+i(f)n-i) 

+a'_T3 exp(27t(/?„_i0„+i) + a'^{n)T^ exp{2'yt(f)n+i(pn-i)} (4.18) 

where 

a± = Vn- 1 ± Vn + 1 (4.19) 

and 

T[ = [a.+l(t)-^^n+l(t)][Cn™l(t)+i.Sn-l(t)] 

•[x„+i(t) - zy;+i(t)][x„_i(i) - zy„_i(i)] (4.20) 

= [C„+i(t)-i5„+i(t)][C„_i(t)-i5„_i(t)] 

•[X,+i(t) - tYn+imXn-i{t) + ^yn-l(^)] (4.21) 

= [C„+i(t)+z5„+i(i)][C„_i(i) + z5„_i(i)] 

•[X„+i(t) + iF„H.i(t)][X„_i(t) - iF„_i(t)] (4.22) 

Ti = [C„+i(t)+25„+i(t)][C„_i(t)-z5„_i(t)] 

•[X„+i(t) + zF„+i(t)][X„_i(t) + iY„_,{t)] (4.23) 

Substituting the explicit expressions of the functions Cn{t), Sn{t), X„(t), and l^(t) 
into the above equations, from Eq.(4.18) we find that 

(a^e^^'"*) ^ jJ2 VnQnQ*n-2 " {a'+nexp[iXa'_{n)t] exp[-26'_(n)7t] 

4 n=0 

+a'^n exp[— iAa'_ {n)t] exp[— 6'^ (^)7^] 
+a'_{n) exp[iXa'_^_{n)t\ exp[— c'_(n)7t] 

+a'_{n) exp[-i\a'^{n)t] exp[-c'_^(n)7t]| (4.24) 
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where 



b'^{n) = [2a;± Ao'_(n)]2, 



4(n) = [2uj ± Xa'^{n)]^ 



(4.25) 



So far we have completed calculations of all expectation values needed in the squeezing 
parameters Eqs.(4.8) and (4.9). It is straightforward to obtain the two squeezing param- 
eters and 5*2 through the simple substitution of the expectation values (4.3), (4.11), 
and (4.24) into Eqs.(4.8) and (4.9). 

We now analyze influence of the dissipation on the quadrature squeezing of the cavity 
field. It is a well-known fact that without dissipation, the field in the JCM governed by 
the von Neumann equation exhibts quardature squeezing [27,28]. Prom Eq.(2.6) we can 
see that when the damping vanishes, it reduce to the coventional von Neumann equation. 
This means that in the JCM with dissipation there exists the quadrature squeezing of the 
field when 7 — >■ 0. Then, taking into account Eqs.(4.3), (4.11) and (4.24), from Eqs.(4.8) 
and (4.9) we can see that the damping term in equation (2.6) leads to the appearance 
of decay factors in each term in the expressions of the squeezing parameters. Thus, each 
term of 5*1 and S2 decays with the decrease of the damping parameter 7. In particular, 
at an arbitrary time t, 7t 1, we find that 



which means that the quadrature squeezing of the field vanishes with the increase of 
the damping. Therefore, we can conclude that the dissipation suppresses the quadrature 
squeezing of the field in the JCM. 

4.3 Sub-Poissonian photon statistics 

A good measure of the extent to which the photon statistics of a state of the light is 
sub-Poissonian is the Mandel's Q parameter [29] defined by 



which characterizes the departure from Poissonian photon statistics. The Mandel param- 
eter vanishes for a Poissonian distribution. When Q < 0, the photons are sub-Poissonian 
while for Q > the photons are super-Poissonian. 



Si — S2 — 



2|z|^ + l > 



(4.26) 



(n^) - {(Any) 
(n) 



(4.27) 
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The expectation value (n) in (4.27) has been evaluated in Eq.(4.3). And the expecta- 
tion value (n^) can be expressed as 

oo 

(n^) = Y,n'^p{n,t) (4.28) 

n=0 

Prom Eqs.(4.2) and(4.28), it follows that 

1 _ °° (2n -\- 1 . 

(n") = (1 + - ^e-" ^ ^ , ^ (n)" exp[-4A'7t(n + 1)] cos[2AiV^H^] (4.29) 

2 n=0 ^• 

where n — 

It is straightforward to get the exact expression of the Mandel's parameter Q through 
the substitution of Eqs.(4.28) and (4.29) into Eq.(4.27). 

As is known, when the damping vanishes, i.e., 7 = 0, the cavity field in the JCM 
can exhibit sub-Poissonian distribtuion [30] in the time evolution governed by the von 
Neumann equation. Moreover, from Eqs.(4.3) and (4.29) we can see that for an arbitrary 
given time t, when X^jt ^ 1, we find that 

(n) = n + i (n^) = (1 + n)^ (4.30) 

which leads to 

«-5^>0 (4.31) 

This indicates that the state of the cavity field is super-Piossonian. Therefore, we can 
conclude that the state of the field in the JCM changes from sub-Poissonian (or Poissonian) 
distribution to the super-Poissonian distribution with the increase of the damping. 

5 Conclusions 

In this paper we have presented an alternative approach to deal with the JCM with the 
dissipation at finite temperature by introducing a new master equation for the reduced 
density operator of the system under the Born-Markovian approximations. We have found 
the analytic solution of the master eqution for the JCM when initially the field is in a 
coherent state and the atom in its excited state. We have studied dynamical properties of 
the atom and the field in the JCM in detail. We have obtained analytic expressions of the 
atomic inversion the atomic dipole moment and the atomic entropy. We have investigated 
infiuence of the dissipation on nonclassical effects in the JCM, such as revivals of the 
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atomic inversion, oscillations of the photon- number distribution, the quadrature sqeezing 
of the cavity field, and sub-Poissonian photon statistics. In particular, we have shown that 
the dissipation suppresses nonclassical effects in the JCM. In our model, the damping 7 
linearly increases with the environmental temperature, therefore these nonclassical effects 
deteriorate with increasing the environmental temperature. 
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